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Abstract 

We consider modified t — J model with minimum of single-hole dispersion at 
the points (0, ±7r), (±7r,0). It is shown that two holes on antiferromagnetic 
background produce a bound state which properties strongly differs from the 
states known in the unmodified t — J model. The bound state is d-wave, it has 
four nodes on the face of the magnetic Brillouin zone. However, in the coordinate 
representation it looks like as usual s-wave. 

PACS: 74.20.-z, 75.30.Ds, 75.50.Ee 



It is widely believed that essential physics of high-T c Cooper Oxide superconductors is 
described by t — J model on two dimensional square lattice. It is well known that at half-filling 
(one hole per site) this model is equivalent to the Heisenberg antiferromagnetic insulator with 
the long-range antiferromagnetic order in the ground state. However any finite doping by holes 
destroys the long-range order (see e.g. Refs. HrH). The problem of two- hole bound state is 
closely related to the problem of superconducting pairing. However it is substantially simpler 
because no destruction of the antiferromagnetic order is to be considered, the two holes can 
not destroy the well denned antiferromagnetic background. 

Single hole (single hole above half- filling) properties in the t-J model are by now well 
established. A hole is an object with complex structure due to virtual admixture of spin 
excitations. The bottom of the hole band is at the momenta k = (±7r/2, ±7r/2). The dispersion 
is almost degenerate along the face of the magnetic Brillouin zone 7k = ^(cosk x + cosk y ) = 
(see e.g. Refs. ^-15 ). So the dispersion at the minima k = (±7r/2, ±7r/2) is very anisotropic: 
the effective mass is very large along the face of magnetic Brillouin zone, in contrast it is equal 
to approximately two electron masses in the perpendicular direction. The mass ratio is about 6 
fi~2yi5[| . A hole can be characterized by the z-projection of spin S z = ±|, where z is a direction 
of the staggered magnetization. 



In the t — J model two holes with opposite S z produce an infinite set of bound states. 
Classification of these states as well as calculation of their binding energies is given in our work 
[IB]. The interesting peculiarity of these states is an unusual relation between the wave function 
in k-representation, parity, and the wave function in coordinate representation. For example 
the state which is usually called d-wave possesses four nodes on the face of magnetic Brillouin 
zone and has positive parity, but in the coordinate representation it looks like a p-wave. This 
peculiarity is due to the existence of the antiferromagnetic background. For t/J < 2 — 3 there 
are also the short-range bound states whose size is close to a lattice spacing. They belong to d- 
and p-wave symmetry according to conventional classification by the number of nodes on the 
face of magnetic Brillouin zone. They are known as well from computer simulations |L7|-|2l] 
and from variational calculations ||22|| . Usually it is convenient to consider these latter states 
in the momentum representation. That is why the unusual properties of the wave function in 
the momentum and coordinate representations are not so interesting for these deep states. 

The pure t — J model gives a Fermi surface with small hole pockets centred at (±7r/2, ±7r/2). 
The mentioned above properties of two-hole bound states are essentially based on the fact that 
minima of single-hole dispersion are at these points. However recent experimental data on 
angle-resolved photoemission (see e.g. Refs. p3"| , |24|| ) demonstrate the existence of large Fermi 
surface with minima at the points (±7r,0), (0, ±7r). Let us incorporate these data into the 
t — J model. With this purpose let us include into the Hamiltonian a small hopping amplitude 
along the plaquette diagonals (compare with Ref. p5|). So let us consider t — t' — J model. 
We will consider very small values of diagonal hopping (£' < O.li), therefore all properties 
of hole-hole and hole spin-wave interactions remain practically unchanged in comparison with 
pure t — J model. However this small t' is enough to shift minima of single-hole dispersion from 
(±7r/2, ±7r/2) to (±7r, 0), (0, ±7r), and we will see that it drastically changes the properties of 
two-hole bound states. In the modified model there is only one exponentially shallow bound 
state of d-wave symmetry having four nodes on the face of the Brillouin zone. In the coordinate 
representation the wave-function of this state possesses no nodes which makes it look like s- 
wave. 

In t — J model for t < 5 the hole dispersion can be well approximated by the following 



analytical expression HpE]] (we set below J 



e k = - v/A 2 /4 + 4t 2 (l +y)- 4t 2 (x + y) 7 £ + -&(cos k x - cos k y ) 2 , 



where A w 1.33, x ~ 0.56, y ~ 0.14. The parameters A,x,y are some combinations of the 



ground state spin correlators [14]]. The fa term gives the dispersion along the face of the 
Brillouin zone. It is rather small and according to Refs. ||i2| , |15| can be estimated as 02 ~ O.lt 
at A/4 < t < 5. The minima of the dispersion ([T]) are at the points (±7r/2, ±7r/2), and the 
mass ratio for different directions at these points is about 6. 

In t — t' — J model we have to add to the dispersion ([!]) the diagonal hopping term 
4t'g/-j cos k x cos ky. The effective hopping parameter t' e ^ differs from the bare one due to the 
dressing of the hole by spin- wave excitations. For t' <C t and A/4 < t < 5 one gets the es- 
timation t' e ff ~ 0.5t'. The correction to dispersion can be represented as ^t' e ^ cos k x cos k y = 
4*e//[7k — |( cos ^ — cos ky) 2 }. The j 2 . term can be neglected in comparison with that in 
Eq. ([!]), and the second term gives the renormalization of dispersion along the face of zone: 
02 02 — 4t' e ff- So the only difference of the modified t — J model from the pure one is a 
different value of 02- Therefore we can consider (3 2 as independent parameter instead of t'. For 
positive fa the dispersion minima are at (±7r/2, ±7r/2), while for negative (5 2 the dispersion 
minima are at the points (±7r, 0),(0, ±7r). 

The state of two holes with zero total momentum and zero projection of spin is of the form: 
\2h) = Skfl , k^ikT^-kj.|0). Here |0) is a wave function of the quantum antiferromagnet (ground 
state of Heisenberg model), is a creation operator of a dressed hole with dispersion (JJ) and 
a — ±|. To avoid misunderstanding note that a number of magnons in a dressed hole is not 
limited. Very roughly at t/J — 3 the weight of the bare hole in is about 25%, the weight 
of configurations "bare hole + 1 magnon" is ~ 50%, and of configurations "bare hole + 2 or 
more magnons" is ~ 25%. The wave-function satisfies the Bethe-Salpeter equation which is 
convenient to write for Xk = (E — 2ek)gk- 

Tkk'Xk' / n \ 



Xk = J2 



k , E — 2e k ' 



Summation over k' is restricted inside the magnetic Brillouin zone (7k' > 0). The kernel Tkk' 



was derived in the work [Jig] . It is the sum of the single spin- wave exchange contribution and 
the "contact" part: Y = T sw + T contact . The single spin- wave exchange part is of the form 



2 (7-k^q + 7k'^q)(7-k'Mq + 7 k Vq) 



r ^ = _ 16/a v ,-^ q ■ ,^ q/v ,-^ ,~w t (3) 
-fc — 6 k ' — €k — Co> q 



where q = k' + k is the momentum of virtual spin- wave, cu q = 2^/1 — 7q is the frequency of 
this wave, u q = + \ and t> q = —sign(jq)^J-^ — \ are the parameters of the Bogoliubov 
transformation diagonalizing spin-wave Hamiltonian. Hole spin-wave coupling constant / is a 



function of t. It was calculated in the work [26[|. The effective interaction (pi) is given by the 



second order of perturbation theory in quasihole-spin-wave interaction. Let us stress that even 
for t > J the quasihole-spin-wave interaction has the same form as for t <C J (i.e. as for bare 
hole operators) with an added renormalization factor which is of the order J / 1 for t^> J. This 
remarkable property of the t-J model is due to the absence of the single-loop correction to the 
hole-spin- wave vertex. It was first stated implicitly in the work ||. In Refs. |T2| , p!3| , |26[| it was 
explicitly demonstrated that the vertex corrections with different kinematic structure are of the 
order of few percent at t/J rs 3. There is also a weak g-dependence of the coupling constant 



/. The value of / which we use in the present work was calculated in the Ref. [26| in the long 
wave length limit: q = 0. At q ~ tt, the factor / is 17-20 % bigger than at q = 0, see the 
discussion between Eqs. (13) and (14) of Ref. [flE}] . This q-dependence does not influence the 
qualitative results of the present work. However it could be important for some quantitative 
results, see conclusion. 

For the states under consideration the energy E in the denominator @ is negligible because 
\E\ <C ek' + ek + c<j q . So we can set in (|3]) E = 0. Then Eq.(||) is reduced to the Schroedinger 
equation. The contact hole- hole interaction is of the form |2l?| , |TE| ] : T contact « 8A7 k _k'- The 
vertices T sw and T contact have been derived for pure t — J model. Here we consider t — t' — J 
model. However t' ~ O.li -C t and therefore the corrections are small. Actually they are very 
small ~ (f /t) 2 because the interactions are due to the spin-flip processes but diagonal hopping 
does not result in the spin-flip. So the interactions T sw and T contact derived for the pure t — J 
model are valid for modified model as well. 

For small t the parameter A in r contac t approaches the value -0.25, which gives the well known 
hole-hole attraction induced by the reduction of the number of missing antiferromagnetic links 
in the problem. However for realistic superconductors t ~ 3 (see e.g. Refs. p?| , p8| ). In this 
region of t the parameter A vanishes due to the hopping contribution |p2| . So this mechanism 
of attraction is switched off. In contrast the spin-wave exchange mechanism (§) is negligible for 
small t where / cx t, and it is the most important at large t where / approaches a t independent 
constant. We are interested mainly in "physical" values of t: t w 3. Therefore in the present 
work we neglect contact interaction (A = 0) and consider only the spin- wave exchange (H). 

Let us concentrate here on the lowest waves of positive parity: d-wave and g-wave. The 
symmetries of corresponding wave functions are shown in Fig.l. It was demonstrated in the 



Ref. [16| that in the pure t — J model (positive ^2) there is an infinite set of bound states in 
each of these waves, and binding energies of the states increase with decrease of fii- Therefore 
consider here first the case /?2 = 0. We found no analytical solution for this case but numerical 



integration of the Eq.flU) is straightforward. The binding energies of four lowest bound states of 
d- and g-wave symmetry at t — 3, / = 1.79, $2 = are presented in Table 1. In the same Table 



we present the values of the radii of bound states in coordinate representation r rms = y(r 2 ), 
where (r 2 ) = / r 2 ip 2 {v)d 2 r = J2k l^f | 2 - From Table 1 we see that at fa — the short-range 
d-wave bound state (r = 1.67) exists at t — 3. Recall that in the pure t — J model it vanishes 
att«2-3 

The negative value of fa changes the situation drastically. At fa < fa g nt the g-wave bound 
state disappears. The critical value is very small. For t — 3, / = 1.79 it is fap 1 = —0.02. 
In the d-wave at fa < f3J lt only one shallow bound state survives. For t = 3, / = 1.79 
the d-wave critical value is faj %t = —0.013. Let us elucidate the structure of this shallow d- 
wave state. Bound state wave function is concentrated around the minima of hole-dispersion 
(0, ±7r), (±7T, 0). Therefore it is convenient to shift the center of the Brillouin zone to one of 
these points. We transfer the usual Brillouin zone centered at k = (0, 0) (solid square at Fig.2) 
to the zone centred at k c = (0, n) (dashed square at Fig.2). This can be done by a combination 
of translations by the vectors of inverse magnetic lattice g = (±7r, ±7r) which are shown in Fig.2. 
The interaction T changes the sign under the translation by g: Tk+g,^ = r^k'+g = — Tkk'- This 
is a consequence of an existence of two magnetic sublattices. One can verify the change of 
sign using explicit formula ([$]) for T sw . The wave function obeys similar relation #k+g = — 9k, 
because it is a solution of Eq.@. Due to this property starting from a d-wave function in the 
usual Brillouin zone (four nodes on the face of zone) after translation (Fig.2) we get constant 
sign wave function inside the shifted zone. The wave function vanishes at the face of this zone. 

For convenience let us define p be a deviation of a momentum from the center of shifted 
zone: p = k — k c , and r p , = r p+ k c ,p'+k c - In the modified model the region of long distances 
corresponds to small momenta p: r ~ p^ 1 ^> 1. Using this fact one can easily prove that 
the Fourier transformation of T pp , gives in coordinate space the short-range potential. This 
situation is quite different from the pure t — J model where the region of long distances corre- 
sponds to small deviation of k from (±7r/2, ±7r/2) and therefore the same interaction gives the 
long-range potential |L6 . 



We consider shallow bound state. Therefore let us solve Eq.(||) using perturbation theory. 
In zero approximation — 1- The integral in the right-hand side of equation (^) is logarith- 
mically divergent at p' = because near this point e p / w const + j\fa\p' 2 . Calculating this 
integral with logarithmic accuracy and taking into account a consistency condition at p = 
one can get the well known formula for exponentially shallow s-wave bound state in two 



dimensions : E ~ — 1/? 2 | • p„ ax • exp(— An^/T^). However in our case Tqq = and therefore 
this approach should be improved. In order to do this let us first evaluate from (fj) a correction 
containing the big logarithm: Xp^ = ~^p o ' A where 

d 2 V 1 ^fMnaA ^ 



2e p -E(2n) 2 4n\f3 2 \ \ E 
Pmax ~ 1- Note that xp=o = 0- Equation @ gives 



x^ + xg^E T l p ' - ■ (5) 



(J.) ,2! = y> ' pP^ 1 ~ V ) 

p' V 



Condition of consistency at small p (i.e. at large distances) xp=o = together with Eq.(|) 
gives the equation for the binding energy 

1 = Q + R-L, (6) 



Q = -E 



1 0,P 

p 2e p : 



d = 1 o.p 1 p,o 
p zt p 

Summation over p is restricted inside the dashed square in Fig.2. From Eq.(||) one finds 

w Ifll 2 ( 4vr^2|(l-g)\ 

£ = -\H -Vmax-^vy ^ y (7) 

So we get exponentially shallow bound state. The integrals Q and R can be easily calculated. 
They are /^-dependent. For t = 3, / = 1.79 and —0.36 < 02 < their typical values are: 
0.3 < Q < 0.7, 0.2 < R < 0.4. The values of exponent in (g) for different fa and t = 3, 
/ = 1.79 are presented in Table 2. In the same Table we present ln^/i?) with E obtained by 
exact numerical integration of Eq.([|). The agreement of analytical formula ([7]) with the exact 
result is perfect. From comparison we find p ma x ~ 1.2. 

We have found exponentially shallow two-hole bound state in modified t — J model having 
minima of the single-hole dispersion at the points (0, ±7r), (±7T, 0). There is no such state in 
the pure t — J model. Thus state behaves as s-wave in coordinate space and also as s-wave in 
momentum space if the rotations around the point k c = (0, 7r) are considered. According to 
the standard classification this state possesses the symmetry of d-wave, - it has four nodes on 
the face of the usual magnetic Brillouin zone centered at k c = (0,0). The size of this state is 



large r rms = yj (r 2 ) = J^^j E\ ^> 1 justifying the analytical formulas used to describe it, but 
numerically formula (0) remains valid even if r rms w 2 when the binding energy is not small. 



The factor R in the exponent of (0) quadratically depends on the interaction potential V. 
This causes extremely strong dependence of the binding energy on the hole spin- wave coupling 
constant /: E ~ exp(— const / f 4 ) . For numerical estimations we have used "infrared" value of 
the coupling constant /. However, it is evident from Eq.@ that the main contribution comes 
from large p, where / is by 17-20% bigger. So numerically we have substantially underestimated 
the binding energy, and account of q dependence of running coupling constant / requires further 
calculations. 

The implication of the found bound state for superconducting pairing will be considered 
elsewhere. However we have to note here that starting from the spin-wave exchange in modified 
t — J model we get the effective short-range attraction. This supports the scenario suggested 



in the recent work [2^] where a short range attraction was introduced by hands. 
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TABLES 



TABLE I. Binding energies (in units of J) and radii of the four lowest bound states in d- and 
g-wave at /3 2 = 0, t = 3, / = 1.79. 





n 


1 


2 


3 


4 


d-wave 


-E 


0.106 


0.17- lO" 1 


0.70 • 10~ 2 


0.39 • 10~ 2 




Trms 


1.67 


5.2 


7.8 


10 


g-wave 


-E 


0.33 • 10- 1 


o.io- io- 1 


0.51 • 10~ 2 


0.31 • 10" 2 




^rms 


4.1 


6.8 


9.3 


12 



TABLE II. Comparison of numerical and analytical calculations for different (t = 3, f = 1.79): 
The value of \n.{(3 2 /E) with E obtained by numerical integration of the Schroedinger equation, and 
the value of exponent in equation (7). 



02 


-0.1 


-0.2 


-0.3 


Hfc/E) 

47r|/3 2 |(l-Q) 
R 


1.0 
1.3 


4.5 
4.3 


10.3 
9.9 



FIGURE CAPTIONS 

FIG. 1. Symmetry of the bound state wave function for (a) d-wave and (b) g-wave. 



FIG. 2. Transition from the usual Brillouin zone 7k > (solid square) to the zone centred 
at k c = (0,7r) (dashed square). 



